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ABSTRACT 

We present the analysis of the Chandra X-ray Observatory observations of the eccentric 7 -ray 
binary PSRB1259-63/LS 2883. The analysis shows that the extended X-ray feature seen in previous 
observations is still moving away from the binary with an average projected velocity of « 0.07c and 
shows a hint of acceleration. The spectrum of the feature appears to be hard (photon index P « 0.8) 
with no sign of softening compared to previously measured values. We interpret it as a clump of plasma 
ejected from the binary through the interaction of the pulsar with the decretion disk of the 0 -star 
around periastron passage. We suggest that the clump is moving in the unshocked relativistic pulsar 
wind (PW), which can accelerate the clump. Its X-ray emission can be interpreted as synchrotron 
radiation of the PW shocked by the collision with the clump. 

Subject headings: pulsars: individual (B1259-63) — X-rays: binaries — stars: neutron — stars: 
early-type — stars: individual (LS 2883) 


1. INTRODUCTION 

High-mass 7 -ray binaries (HMGBs) consist of a com¬ 
pact object (i.e., a black hole or a neutron star) and 
a massive (early-B or late-0 type) star. Only one of 
the few known HMGBs, B1259-63/LS 2883 (hereafter 
B1259), has a compact object that has been detected as a 
radio pulsar. The pulsar has the following properties: pe¬ 
riod P=47.8 ms, characteristic age r = P/2P = 330 kyr, 
spin-down power if=8.3xl0^® erg s“*^, r nagnetic field B 
= 3.3 X 10*^*^ G, and distance d « 2.3 kpc (Johnston et al. 


19921. The pulsar’s companion is a fast-rotating O star, 
with mass » 30M(t, an d luminosity L* = 6.3 x 
( [Negueruela et ah 2011), and an equatorial decretion 
dislQ which is inclined at an angle of » 35° to the or - 
bital plane (Melatos et al. 1995| Shannon et al. 2014). 
The system has an orbital period i^rb = 1236.7 days, an 
eccentricity e = 0.87, a semi- major axis a » 7AU, an d 
an inclination angle i « 23° (Negueruela et al. 2011). 


The pu lsar is eclipsed for abo ut one month around peri- 
astron ( Johnston et al.|[2005 ) when it is moving through 
the decretion disk and dense wind of its companion star. 
During the pulsar eclipse, peaks of nonpulsed radio emis¬ 
sion were observed, centered at « 5 days before perias¬ 
tron and « 20 days after periastron, with the rel ative 
strengths of the peaks varying from cycle to cycle. |Ball| 
et al. (1999) interpreted it is as synchrotron radiation of 
electrons accelerated during the pulsar passages through 
the decretion disk. 

B1259 has been extensively studied in X-rays and 7 - 
rays. The X-ray emission from B1259 showed no pulsa- 


^ Department of Astronomy & Astrophysics, Pennsylvania 
State University, 525 Davey Lab, University Park, PA 16802, 
USA; ggpl@psu.edu 

^ George Washington University, 105 Corcoran Hall, Washing¬ 
ton, DC 20052, USA 
^ Toronto, ON, Canada 

^ Decretion disks (also known as excretion disks) are viscous 
gaseous disks around rapidly rotating O or B stars (Be stars) in 
which stellar matter slowly moves outward and rotates with nearly 
Keplerian velocity in the equatorial plane (see Rivinius et al. 2013 
for a recent review). 


tions. It was interpreted as synchrotron radiation from 
relativistic pulsar wind (PW) leptons accelerated at the 


shock between the PW and the massive star outflow (Ta- 


vani & Arons 1997). Multi-epoch measurements with 


numerous X-ray observatories revealed binary phase de¬ 
pendences of the flux, photon index, and absor ption col¬ 
umn density (Ghernyakova et al. 2006 2009). Within 
the same orbi'tal cycle, the strongest variations are seen 
around periastron. In particular, the X-ray light curve 
shows two peaks corresponding to first and second pas¬ 
sages of the pulsar through the disk (« 20 — 4 days before 
periastron and « 10—50 days after periastron), while the 
hydrogen column density, « 3 x 10^*^ cm“^ during 
most of the cycle, increases by a factor of 2 during and 
between the two disk passages, conhrming the presence 
of additional absorbing matter at these binary phases. 

This unique system is one of the few HMGBs whose 
flaring in the GeV energy range was seen by the Fermi 
Large Area Telescope. The GeV flares of B1259 were 
detected soon after the pulsar passed through periastron 
in 2010 and 2014, while no GeV emission was seen when 


2011 Tam et al. 


it was far f rom periastron (Abdo et al. 

20111 2015). The flares started « 20 — ll5 days after the 
epochs of periastron and lasted roughly 1-2 months. 
This GeV emission could be synchrotron radiation from 
PW electrons accelerated to energies ^ 100 TeV in the 
intra-binary shock, or inverse Gompton (IG) scattering of 
photons emitted by the mass ive companion off electrons 
with energies ^1 — 10 GeV (lAbdo et al.||2011|). 

In the TeV energy range, B1259 has been obse 


"observed by 


the H.E.S.S. observatory (Aharonian et al. 2005 2009). 
Variable TeV emission was detected around the perias¬ 
tron passages of 2004 and 2007, with light curves similar 
to those of the nonpulsed radio emission. No signatures 
of TeV emission enhancement were detected at the times 
of GeV flares after the 2010 and 2014 pe riastron passages 
(Abdoetal. 2011 Malyshev et al. 2014). Several models 
have been developed to explain the TeV emission, such as 
IC sca ttering of stellar photons off ultra-relativistic elec- 


trons (Aharonian et al. 2009 
and hadronic interaction of t 


Khangulyan et al. 2007) 
le PW with the dense de- 
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cretion disk (Aharonian et al.||2009). 

The size oi the B125y binary is too small to resolve it 
in X-rays, but the high angular resolution of the Chan¬ 
dra X-ray Observatory (CXO) allows one to look for ex¬ 
tended emission from the pulsar wind nebula (PWN) or 
from matter that could be expelled from the binary. First 
detection of exte nded emission assoc iated with B1259 
was reported by Pavlov et al. (2011). They observed 
B1259 with CXO on 2009 May TT^bsID 1089; 26 ks 
exposure), when the pulsar just had passed the apas- 
tron of 2009 April 5 (see Figure 1), and detected faint, 
asymmetric extended emission seen up to « 4" south- 
southwest from the binary position. This emission was 
interpreted as a PWN blown out of the binary by the 
massive companion’s wind. 



Fig. 1.— Orbital locations of the pulsar at the times of our 
Chandra ACIS observations, marked by their ObsIDs. Our ‘new 
data’ consists of two exposures (ObsIDs 16563 and 16583) sepa¬ 
rated by about 10 hours. Gray areas show the parts of the orbit 
where the pulsar passes through the companion’s equatorial disk. 


Kargaltsev et al.| (2014 referred to as K-l-14 hereafter) 
reported results of two deeper CXO observations (Ob¬ 
sIDs 14205 and 14206; 56 ks exposures). These obser¬ 
vations were carried out on 2011 December 17 and 2013 
May 19, i.e., 370 and 886 days after the periastron of 
2010 December 14. They revealed a variable, extended 
(about 4", ^ 1000 times the binary orbit size) structure, 
which appeared to shift by 1'.'8±0'.'5 between the two ob¬ 
servations. K-l-14 discussed two possible interpretations 
of the observed variability: (1) an extrabinary termina¬ 
tion shock in the PW outflow, with a variable stand-off 
distance, or (2) a clump ejected from the binary and 
moving with a projected velocity vj_ = (0.046 ± 0.013)c, 
at d = 2.3 kpc. The analysis of those data favored the 
scenario in which most of the observed extended emission 
would come from the PW launched from the binary near 
the apastron of 2009. To gain a better understanding of 


the observed phenomenon and test our interpretations, 
we proposed a Director Discretionary Time (DDT) CXO 
observation. Here we report the results of the new ob¬ 
servation, analyze it jointly with the two previous ob¬ 
servations (taken within the same orbital cycle between 
the periastrons of 2010 and 2014), and suggest a new 
interpretation of the extended structure. 

2. OBSERVATIONS AND DATA REDUCTION 

B1259 was observed with the Advanced CCD Imaging 
Spectrometer (ACIS) on-board CXO on 2014 February 
8 and 9 (ObsIDs 16563 and 16583, respectively), with 
a total live time of 57.6 ks. The corresponding orbital 
location of the pulsar (« 1150 days after the peristron of 
2010; 9 = 221°) is shown in Figure 1. The observation 
setup was similar to our previous observations that have 
been analyzed in K-l-14. The target was imaged on the 
front-illuminated ACIS-I3 chip in timed exposure mode, 
the data were telemetered in ‘very faint’ format. A 1/8 
subarray was used to reduce the frame time to 0.4 s and 
mitigate the effect of pile-up. The highest count rate 
observed from the binary (0.043 counts per frame) cor¬ 
responds to a negligibly small pile-up fraction of < 2%. 

We used the pipeline-produced Level 2 event files for 
the analysis. No episodes of anomalously high back¬ 
ground rates occurred in the observations. To minimize 
the background contribution, we filtered the events by 
limiting the photon energies to 0.5 — 8 keV. The detector 
response for spectral analysis was produced with CIAO 
(ver. 4.6) tools, following the standard procedure and us¬ 
ing the calibration database CALDB 4.5.9. The spectral 
fitting was done with XSPEC (ver. 12.8.0). 

3. DATA ANALYSIS 

3.1. Apparent motion of the extended feature 

Since the two new observations were carried out less 
than one day apart, we merged the data for purposes 
of image analysis. The merged image clearly reveals a 
‘blob’ of extended emission with a very faint westward 
extension. The appearance of the extended feature has 
changed substantially between the precediM observation 
ObsID 14206 and the new data (see Figure!^. Unlike the 
earlier images, in which the feature resembled an arc, the 
feature’s shape is more round in the new image. It is 
possible that the fainter part of the feature had simply 
become too faint to be detected as a result of a decrease 
in the surface brightness. 

The new data allow us to further track the motion of 
the extended feature reported in K-l-14. Using the CIAO 
celldetect tool, we measured the positions of the bi¬ 
nary and the extended feature together with the corre¬ 
sponding la uncertaintieFl The 3'.'1 ± Of'7 shift of the 
blob, with respect to the position of the arc-like struc¬ 
ture in the ObsID 14206 image, corresponds to a proper 
motion p = 4'.'2 ± 0'.'9 yr“^ and projected velocity t>_L = 
(0.15±0.04)c, at d = 2.3 kpc. The temporal dependence 
of the radial separation of the extended emission from 

® For the elongated arc-like structure in the earlier images, 
K-l-14 used a different method (see Figure 3 in K-l-14), but running 
celldetect on ObsID 14206 gave consistent results. For ObsID 
14205 celldetect fails to pick out the azimuthally extended struc¬ 
ture because it was too close to the bright binary, so we used the 
peak in the radial count distribution for the blob position and the 
FWHM of the peak for its uncertainty. 
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the binary, at a position angle of 215° (north through 
east), is plotted in FigureThe plot does not show any 
evidence of deceleration, rather it might suggest that the 
blob is moving with an acceleration u_l ^ Aw_L/At ^ 
(0.010 ± 0.04)c/(390d) = 90 ± 40 cms“^. Given the 
low significance of this result, we assumed a constant 
velocity and estimated it by fitting the radial separa¬ 
tions with a straight line, r{t) = fi(t — Iq) -|- tq, where 
to = 801 days is the reference time (all times are counted 
from the periastron of 2010 December 14). The fit yields 
ro = 4"6±(y.'4, fj. = (y.'0056±0'.'0009d-i = 2'.'0±0"3yr-h 
This proper motion corresponds to the projected velocity 
f_L = (0.07 ± 0.01)c at d = 2.3 kpc. The launch time, 
defined as r(tiaunch) — 0, is tiaunch — 21_^gr^ days. This 

range of tiaunch encompasses the periastron and the Fermi 
flare dates. 


3.2. Spectral analysis 

We extracted the spectra from each of the two new ob¬ 
servations and fitted them simultaneously. The extrac¬ 
tion regions are shown in Figure We fit the spectra 
of the compact core (6257 counts in the Oi'O radius aper¬ 
ture) and the blob (58 counts in the 3.6 arcsec^ area) with 
the absorbed power-law model, using the XSPEC phabs 
model for absorption. We obtained the following fitting 
parameters for the core: photon index F = 1.5 ± 0.4, hy¬ 
drogen column density Nh — (4.2 ± 0.3) x 10^^ cm“^, 
and unabsorbed flux = (1.95 ± 0.04) x 10“^^ 

erg cm“^ s“^ (see Table 1 for more details). 

Because of the small number of counts in the blob, we 
used C-statistic[^ (without binning) and fixed iVn to fit its 
spectrum. The fitting parameters are provided in Table 
1. Figured shows the 90% and 99% confidence contours 
in the T-W plane {Af is the power-law normalization) 
for Nh = 3 X 10^^ cm“^. The latter value is smaller 
than the one we obtained from fitting the core spectrum; 
we chose it as the average of the Ah values obtained 
in the fits of the core spectrum in ObsIDs 14205 and 
14206 (K-l-14) in order to exclude additional absorption 
from within the binary itself. Figure H also shows the 
T-Af confidence contours for the extended emission from 
two previous observations, which we re-fit using the C- 
statistic for consistency. We see that the photon index for 
the blob in the new data, F = 0.8 ± 0.4, does not exceed 
those found in the previous observations, F = 1.2 ± 0.2 
and 1.3 ± 0.2 for ObsIDs 14205 and 14206, respectively 
(i.e., the spectral evolution does not show any cooling of 
the blob’s matter). We see from Figure]^ that the flux of 
extended emission has been steadily decreasing. Fitting 
it with an exponential decay law yields a characteristic 
decay time of 540 ± 100 days. 


4. DISCUSSION 

K-l-14 have discussed possible scenarios regarding the 
origin and nature of the observed variable extended emis¬ 
sion. One possibility considered by K-l-14 was the syn¬ 
chrotron emission from the PW that escapes from the 
binary in the direction away from the high-mass star and 
shocks in the extrabin ary medium, similar to PW Ne of 
young solitary pulsars (Kargaltsev fc Pavlov||2008 ). The 


® https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/ XSap- 
pendixStatistics.html 


estimated stand-off distance to such a shock is compara¬ 
ble to the observed separations for an ambient pressure 
Pamb 10“^° dyn cm“^. Because of the large eccentric¬ 
ity of the binary, the pulsar spends most of orbital cy¬ 
cle around apastron, so most of the shocked PW matter 
concentrates around the extension of the binary’s major 
axis in the apastron direction. The separation between 
the extrabinary shock and the pulsar could vary due to 
changes of the ratio of the PW pressure to the ambient 
pressure. However, this interpretation looks somewhat 
artificial now because the new data strongly suggest that 
the separation is steadily increasing, which would require 
a steadily decreasing ambient pressure (or increasing PW 
pressure) on a time scale of ~ 800 days. In addition, the 
required ambient pressure is unrealistically high for a cir- 
cumbinary medium, which is likely filled with a freely 
flowing stellar wind of the massive companion, and the 
extended structure does not have an arc-like appearance 
anymore, hence losing its resemblance to a termination 
shock. 

Another scenario considered by K-l-14 assumed that 
the moving structure is a clump of matter ejected from 
the binary by the pulsar interaction with the decretion 
disk. The new data seem to support this scenario because 
the average projected velocity implies that the structure 
could indeed be launched near periastron. However, as 
K-l-14 pointed out, the fast, steady motion over such a 
large period of time poses a number of problems. If the 
ejected clump consists solely of synchrotron-emitting rel¬ 
ativistic electrons or positrons, possibly confined by the 
magnetic field, it should decelerate very quickly by the 
drag force (assuming it moves through a much slower 
expanding medium, such as a bubble created by the 
companion’s wind). On the other hand, if the clump 
consists solely of electron-ion plasma (e.g., a fraction 
of the decretion disk heated and expelled from the bi¬ 
nary by the interaction with the PW)[^ then its mass 
would be rrici ~ 10 ^® g, much larger th an the plausible 
ma ss of t he disk, mdigk ~ 10^"* — 10^® g ([Chernyakova et 


al. 


2014 


and references therein), or the mass MPorb = 
6.7 X 10^®(M/10“® M 0 /yr) g lost by the high-mass com¬ 
panion during one orbital cycle. In addition, its kinetic 
energy, Ed ~ 10 ^^ erg, would exceed the energy that 
could be supplied by any conceivable source. 

A mixed scenario, in which the emission is produced by 
relativistic electrons supplied by the PW but the clump’s 
mass is determined by protons/ions from the stellar out¬ 
flow, may be more plausible. However, even in this sce¬ 
nario, the clump’s mass should be a substantial fraction 
of m^isk to overcome the drag force, even at extremely 
low ambient densities (K-l-14). Moreover, it remains un¬ 
clear what is the source of the corresponding clump’s 
kinetic energy, Ed ^ 5 x 10 "‘^(toci/ 10^® g)(u/0.1c)^ erg, 
let alone the acceleration mechanism. If the clump’s 
energy is provided by the loss of a small fraction of 
the pulsar’s spin energy, Egp = 8.6 x 10 '^® erg (for the 
moment of inertia I = 10"^® g cm^), during the disk 
passage, the decrease of the spin energy, AAgp ^ Ad, 
would lead to a perceptible increase of the spin period. 


^ In this case, the clump emission would be thermal 
bremsstrahlung. As shown by K+14, the spectra of the extended 
feature can be fitted by a thermal bremsstrahlung model with 
kT > 10 keV, rie ~ rip ~ 10^ cm“^. 
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Fig. 2. — ACIS images from our previous two observations, OblDs 14205 (top left) and 14206 (top right), and the new data (bottom left). 
The blob of extended emission is clearly traveling away from the binary. The bottom right panel shows the false color image (different colors 
correspond to different observation epochs) demonstrating the motion and changes in morphology. The dashed line connects the binary 
and the centroid of the feature in the new data image. The angular separation between the extended feature and the binary as a function 
of time is shown in Figure]^ The spectral extraction regions (core and extended emission are shown in black and cyan, respectively) on 
new data image. North is up, east is to the left. The radius of the black circle is 0^9; the background extraction region (not shown) 
is an annulus with radii 8^55 ^ r < 19^18. The white crosses mark the positions of the extended emission as obtained with the CIAO 
celldetect tool. The apparent asymmetry of the bight core seen in the merged image of ObsIDs 16563+16583 can, at least partly, be 
attributed to the known Chandra mirror asymmetry. 


AP ^ (P^/47r^J)i?ci ~ 6 X 10“®(mci/10^® g) s. However, 
no such period jumps wer e noticed in the most detailed 
timing analysis available (Shannon et al.||2014). These 
authors mention systematic variations in timing residu¬ 
als close to the 6 periastron passages covered by their 
analysis, but they attribute those to uncorrected varia¬ 
tions in dispersion measure or scattering of the pulsar 
radiation through the disk. Thus, it seems that the spin- 
down due to the pulsar-disk interaction does not look as 
a plausible source of the clump energy. In principle, one 
might assume that a small fraction of pulsar’s kinetic 
energy, Pkin 10'^® erg (near periastron), is somehow 
converted into the clump’s energy in the course of pulsar- 
disk interaction, but a decrease of Pkin should result in 
a decrease of binary period Pb, which is not consistent 


with the results of timing analysi^ Thus, we conclude 
that even the mixed clump scenario looks problematic 
because ejection of a large mass is required to overcome 
the drag force, and it remains unclear what is the energy 
source and how such a mass could be accelerated to a 
velocity ~ 0.1c during the short time of the pulsar-disk 
interaction. 

The drag deceleration and energy deficit problems can 
be alleviated if, instead of the companion’s wind bubble, 
the clump is moving in the unshocked PW whose veloc¬ 
ity, close to the speed of light, significantly exceeds the 


® On the contrary, [Shannon et al.]j2014| found a hint of increase 
of the binary period, with a derivative = (1.4 it 0.7) X 10“® s 
s“^, corresponding to APj Ri 1.5 s per binary cycle, which they 
attribute to mass loss due to the polar wind of LS 2883. 
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Days since periastron passage 


Fig. 3. — Separation of the extended feature from the unresolved 
point source (in arcseconds) as a function of time since the peri¬ 
astron passage on 2010 December 14. The line of best fit and Icr 
upper and lower bounds are also shown. The small red line at 
the bottom shows the time pe riod. 2011 January 14 - March 4, in 
which the GeV flare occurred l |Abdo et al.|20li] |. 


clump’s velocity. This scenario becomes more plausi- 
ble at larger values of the parameter 77 = E/{Mvwc) = 
4A MZgV~\, where M-g is the mass-loss rate in the O- 
star polar wincQ in units of 10 “® Mq yr~^, and is 
the wind’s velocity in units of 1000 km s“^. At 77 > 1, 
the companion’s wind is confined by the PW into a cone 


with half-opening angle a « 30°(4 — (Eich- 

ler &: Usov|l993 ), while the unshocked PW fills the rest 
of circumbinary volume^ Since the pressure in the 
freely flowing (unshocked) polar wind of the high-mass 
companion is proportional to similar to the PW ram 
pressure, Ppw = E/( 47 rcr^) = 2.2 x 10 “^°ry dyn cm“^ 
(where ryj = r/10^^ cm), the PW cannot shock until 
it reaches the termination shock of the stellar wind at 
r = rts, where the wind bubble pressure j umps up by 


a few orders of magnitude (see Figure 1 in Dwarkadas 


20051. This means that the PW shocks at large dis- 
tances from the binary (e.g, at ~ 3' for rts ~ 2 pc), 
where the surface brightness of the shocked PW emis¬ 
sion (PWN) may be too low to be detected. In this 
scenario, the observed extended X-ray emission can be 
due to synchrotron radiation of the PW shocked by the 
collision with the clump of ejected material. The X-ray 
luminosity can be estimated as Tx.ci = VxE{rci/2r)^, 
where r is the distance from the pulsar, r^ is the effec¬ 
tive radius of the clump, and rjx is the X-ray efficiency. 


® The reference value of mass-loss rate, M = 10 ® Mq yr 
is lower than usually adopted for LS 2883, but such a low M is 
consistent with observations of OV stars with luminosities L* < 
1 O®'^L 0 , which show mass losses lower than theoretically predicted 
or expected from extrapolations from more luminous O stars (the 
so-called “weak wind problem”: see Section 5 in the review bvIPulsI 

let al.||2008l) '-' 

A large value of 77 is needed in our scenario because in the op¬ 
posite case the unshocked PW would occupy only a narrow channel 
confined by the stellar wind and rotating around the massive star 
together with the pulsar. Even if a clump of stellar matter gets 
into that channel, it will go out from the channel in a fraction of 
the binary period and will be decelerated by the stellar wind. We 
note that larger rj also imply the intra-binary shock is closer to the 
massive star, which helps to explain the hig h-enerev emission by 
the Compton up-s cattering of stellar photons (| Sierpowska-B artosikj 
|fc Bednarek|2008|l. 


For instance, Tx,ci ~ O.OlrjxE ~ 8 X 10^^77x erg s ^ 
for Tci/r ^ 0.2 for the latest observation. This esti¬ 
mate is consistent with the observed X-ray luminosity, 
Ax,cl = (1-3 ± 0.3) X SI'S at a reasonable 

value for X-ray efficiency, rjx ~ 1.5 x 10“^. This inter¬ 
pretation of the X-ray emission is also consistent with 
the lack of spectral softening (Section 3.2) because cool¬ 
ing of relativistic electrons is compensated by the energy 
supplied by the PW. 

The interaction of the unshocked PW with the ejected 
fragment is not only responsible for its X-ray emis¬ 
sion, but it can also accelerate the fragment. A very 
crude estimate of acceleration is 7 ) Ppv,Am~^^ 
440rj"7^^A(wci/10^^ g)“^ cm s“^, where Ppw is the PW 
ram pressure, A 2 x cm^ is the clump’s cross 

section, and < 1 is the filling factoip^ Interestingly, 
this estimate is consistent with the apparent accelera¬ 
tion estimated (with low significance) in Section 3.1 at 
TTici lO^^^A g (for r = 2 X 10^^ cm, i.e., ^ 7" from the 
pulsar). Ejection of such a small mass should not make 
a measurable effect on pulsar timing. Moreover, the pos¬ 
sibility of gradual acceleration of the ejected clump alle¬ 
viates the problem of a huge acceleration of the clump in 
the course of its ejection caused by pulsar-disk interac¬ 
tion. For instance, for mc\ ^ 10^^ g, the clump’s kinetic 
energy, Ed ^ 4.5 x 10 ^®( 7 ;/ 0 . 1 c)^ erg, is a tiny fraction, 
^ 5 X 10“^, of the pulsar’s spin-down energy loss dur¬ 
ing, e.g., the time elapsed from the 2010 periastron. We 
should note that if the acceleration during our observa¬ 
tions was indeed so large, then the conclusion that the 
clump was launched around the 2010 periastron, derived 
under the assumption of constant speed, becomes ques¬ 
tionable (unless the acceleration was much lower when 
the clump was closer to the binary, e.g., due to a larger 
accelerated mass). The matter could be ejected by the 
pulsar-disk interaction close to an earlier periastron, but 
it is hardly possible to estimate the launch time without 
detailed modeling of the PW-clump interaction, account¬ 
ing for varying acceleration. 

A qualitative scenario for cloud formation and accel¬ 
eration might look as follows. Before entering the disk, 
the pulsar moves in the fast, radially directed compan¬ 
ion’s polar wind, the apex of the intrabinary bow shock 
is on the line connecting the pulsar and the compan¬ 
ion, and the shock position and shape are determined 
by the parameter 77 (Dubus||2013|). When the pulsar en¬ 
ters the disk, the structures of both the disk and the 
intrabinary shock change dramatically. Since the mat¬ 
ter in the disk moves in almost circular orbits around 


the compan ion with nearly Keplerian velocity (Rivinius 
et al.||2013), which is comparable to the orbital velocity 
of the pulsar, the relative velocity of the pulsar and the 
disk matter, Vj-ei = Vdisk — Vpsr, is not directed radially, 
and the shock apex is not on the pulsar-compani on line 
(see an illustration in jKhangulyan et al.| |2012|). The 
(minimum) distance of the shock from the pulsar can 
be estimated as Vg ~ i?^/^[47rcpdisk(R)]~^^^|vrei(R)|“^, 
where pdisk{R) ~ po{R*/R)^ is the density in the mid- 


The ‘clump’ can consist of a few smaller fragments, which we 
cannot resolve in the image. These fragments could be launched 
from the binary in slightly different directions, which could explain 
the large size of the observed extended feature. 
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TABLE 1 

Spectral fit parameters for the core and extended emission in three CXO ACIS-I observations 


ObsID 

MJD 

ga 

deg 

AA 

days 

Exp.^ 

ks 

els'* 

Hjbs® 
10“^^ cgs 

TP f 

^corr 

10“^^ cgs 

Nh 

lO^i cm-2 

r 

10-"* 

(arcsec^) 

X^/dof 

10089 

54965 

182 

667 

25.6 

1825 

139(5) 

158(6) 

1.5(7) 

1.51(10) 

2.3(3) 

2.5 

20.9/30 






61 

2.8(1) 

3.1(1) 

1.5* 

1.3(5) 

0.039(16) 

22.1 

3.1/9 

14205 

55912 

169 

370 

56.3 

6551 

249(4) 

296(5) 

2.9(3) 

1.39(5) 

3.8(2) 

2.5 

76.6/87 






343 

8.5(5) 

9.2(7) 

2.9* 

1.2(1) 

0.10(1) 

22.1 

36.24%k 

14206 

56431 

192 

886 

56.3 

4162 

137(5) 

176(7) 

3.1(3) 

1.68(6) 

3.07(2) 

2.5 

146/169 






144 

3.6(4) 

3.9(6) 

3.1* 

1.3(2) 

0.052(14) 

12.8 

62.96%!^ 

New^ 

56696 

221 

1151 

57.6 

6257 

149(5) 

195(4) 

4.2(3) 

1.5(4) 

3.0(1) 

2.5 

242.5/204 






58 

1.9(4) 

2.0(4) 

3.0*-i 

0.8(4) 

0.013(7) 

3.6 

54.12%'^ 


Note. — For each ObsID the upper and the lower rows correspond to the partially resolved core and the extended emission, 
respectively.The Icr uncertainties are shown in parentheses (see also Figure 4). The XSPEC extinction model phabs was used throughout. 
Fluxes and counts are in the 0.5—8keV range, corrected for finite aperture size for the core emission. An asterisk indicates that the 
extinction column was fixed to the value of the corresponding point source fit. 

^ True anomaly counted from periastron. 

^ Days since latest preceded periastron. 

^ Exposure corrected for deadtime. 

^ Total (gross) counts. 

® Observed flux. 

^ Extinction corrected flux. 

s Normalization in photonss~^ cm“^ keV“^ at 1 keV. 

^ Area of the extraction region. 

^ Obsids 16563 and 16583 fit simultaneously. 

j Fit was done, for extended emission, using the fixed value obtained from fitting the core around apastron passage to exclude any 
absorption from within the binary. 

^ C-statistics were used for fitting the extended emisssion spectra across all observations for consistency. This percentage is the Null 
Hypothesis probability 1 — P, where V is the probability given by the XSPEC goodness command with 10,000 runs. 



Fig. 4.— Confidence contours (90%, and 99%) in the P-A/” plane 
for the ACIS spectra of the extended nebula (see Figure and 
K+14 for region definitions) computed for the absorbed PL model. 
The PL normalization J\f is in units of 10“^ photons cm“^ s~^ 
keV~^ at 1 keV. The dashed lines are lines of constant unabsorbed 
flux in the 0.5—8 keV band (labeled in ergcm”^ s~^). The contours 
for observations 14205, 14206 and the new data were fit with Nh 
fixed at the values of A^h = 2.9, 3.1, 3.0 x 10^^ cm“^, respectively. 


plane of the disk in the pulsar vicinity (at the distance 
R from the companion), ^ IORq is the compan- 
ion’s radius. Using the same disk model as |Chernyakova| 
eTH] dMIt , we obtain the shock radius ^ 1.5 x 


10Tpdisk/10“^®gcm“^)“^/^(?;i.ei/100kms“^)“^ cm, for 
the first crossing of the disk (i? ^ 40i?* ~ 2AU). Since 
Cg exceeds (or is comparable to) the vertical size of the 
disk, H{R = 40i?*) ~ 5 x 10^^ cm, the disk should be 
significantly disrupted by the first pulsar passage, which 


is consistent with nume rical simulations (Okazaki et al. 
2011 Takata et al.|2012). In the second disk passage fur¬ 


ther fragmentation ot the disk can occur. The collision 
of the PW with disk matter fragments would result in 
shocks and y-ray flares from the shocked PW, as well as 
entrainment of the clumps in the PW flow, which would 
blow them out of the binary and eventually accelerate up 
to the observed velocity ~ 0.1c. The details of the mass 
loading and ejection processes are likely very complex 
and should be investigated with numerical simulations. 

A possible mechanism of matter ejection from the bi¬ 
nary could be connected with episodic accretion of the 
disk matter onto the neutron star magnetospher^pj If 
there were no PW and associated shock, the neutron 
star moving in the decretion disk could gravitation¬ 
ally capture the matter within the cylinder of radius 
i?grav = 2GMp,,v-^ = 3.7 X 10i2(r;,,i/100kms7i)-2 
cm (for Mpsr = I. 4 M 0 ), with the maximum 

possible capture rate mmax = 7ri?gravPdisk?^rei = 
4.4 X 10i6(pdisk/10-i®gcm-3)(z;,,i/100kms-i)-3 g s'k 
However, if the radius Vg of the PW shock exceeds 
i?<r 


(as for the parameters assumed above for the 
shock radius estimate), the matter is deflected by the 
shock out of the capture cylinder. Nevertheles, the cap¬ 
ture becomes possible if i?grav > i-e., if Pdisk > 

1.6 X 10“^®('(;i.ei/100kms“^)^ g cm“^. If the disk den¬ 
sity is indeed so high in the vicinity o f the crossing 


pulsa r (which is not ruled out — see Takata et al. 
2012[), then the ram pressure of the gravitationally at- 


tracted matter {pvg oc r ^/^) overcomes the ram pres¬ 
sure of the PW (ppw cc r“^), and the disk matter 


Accretion onto the neutron star surface would be inconsistent 
with the relatively low X-ray luminosity of B1259 even around 
periastron. 
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streams toward the pulsar with a nearly free-fall ve¬ 
locity, Uff = (2GMpsr/r)^/^, until it reaches the pul¬ 
sar magnetosphere at r ~ Ric = cPj{2TT) — 2.3 x 10® 
cm. At radii smaller than the light cylinder radius Ric 
the accreting matter is stopped by the magnetic field 
pressure oc r“®) in the pulsar magnetosphere, 

at the Alfven radius Ra ~ ~ 


different directions, we have to assume that the radio and 
X-ray features are intrinsically different (e.g., comprised 
of different matter) or that the ejected fragment could 
change the direction of its motion because of the varying 
direction of the accelerating force close to the binary. 

5. CONCLUSIONS 


1.4 X cm, where fi = Bns-Rns ^ 

G cm® is the pulsar’s magnetic moment, and rhn = 
rh/(10^^gs“®). As long as Ra exceeds the corotation 
radius i?cor = {GMP’^jATry^^ = 2.2 x 10^ cm (i.e., at 
any reasonable accretion rate rh < 6 x ipiQ g s~®), accre¬ 


tion occur s in the ‘propeller regime’ (Illarionov & Sun- 
yaev||1975 ), i.e., the accreting matter does not reach the 
neutron star surface. Most of this matter is expelled 
in a disk-shaped outfl ow in the equator i al pla ne by the 


centrifugal forces (see Romanova et al. 2003| for MHD 
simulations) while up to ~ 107o can leave the system 
in a high-velocity jet along the neutron star spin axi^^ 


( Lovelace et al.|[2014[ and references therein), with a ve¬ 
locity up to {GM/RaY^"^ ^ O.Odrhiyc. An upper limit 
on the total mass ejected in a jet during the disk passage 
can be as large as ~ O.lmtdisk 10^^mi7(tdisk/15d) g. 
If the ejected matter gets into the region of unshocked 
PW, it can be further accelerated to the observed veloc¬ 
ity of ~ 0.1c. Thus, accretion of decretion disk matter 
in the propeller regime could, in principle, be responsible 
for ejection of the observed clump, but it would require a 
rather high accretion rate (i.e., a very large disk density). 

Overall, the interpretation of the moving extended X- 
ray feature as a fragment of the decretion disk ejected 
from the binary by the pulsar-disk interaction and mov¬ 
ing in the unshocked PW looks quite plausible. However, 
we should mention a possible contradiction between such 
an interpretation and the detections of extended radio 


Our new CXO observation of the extended object in 
the vicinity of the B1259 binary has shown that the ob¬ 
ject keeps moving away from the binary. The comparison 
with two previous observations shows that its luminosity 
has been steadily decreasing, and its average projected 
velocity was u_l « 0.07c during 780 days covered by the 
observations. If that velocity has remained constant in 
the course of motion, the object was ejected from the 
binary near the 2010 periastron, perhaps by the collision 
of the pulsar with the decretion disk of the high-mass 
companion. If this clump of matter were moving in a 
slowly expanding circumbinary medium, its mass and ki¬ 
netic energy would be uncomfortably large to overcome 
the drag deceleration. Therefore, we suggest the clump 
(likely a fragment of the decretion disk) is moving in the 
unshocked relativistic PW. Its X-ray emission can be in¬ 
terpreted as synchrotron radiation of the PW shocked by 
an interaction with the fragment. The ram pressure of 
the PW can be responsible for clump acceleration, a hint 
of which is seen in our observations. No such phenom¬ 
ena have been observed so far. Further monitoring with 
CXO is needed to confirm or refute this interpretation. 


Support for this work was provided by the National 
Aeronautics and Space Administration through Chan¬ 
dra Awards DD3-14070 and GO2-13085 issued by the 
Chandra X-ray Observatory Genter, which is operated 
by the Smithsonian Astrophysical Observatory for and 


astron ( 

Ghernyakova et al. 2014 

and 21 days after the tion under contract NAS8-03060. We are very grateful 

2UU7 periastron ( 

Moldon et al. 

2011 

). 'ihat emission, to Harvey Tananbaum tor allocating the UUT time to 


and a larger extension toward northwest, was interpreted 
as a cometary tail extending behind the pulsar in the 
northwest direction, quite different from the southwest 
direction of the X-ray feature motion. To reconcile the 


Khangulyan for useful discussions and the anonymous 
referee for careful reading of the paper and useful com¬ 
ments. 
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